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THE MOTION OF A POINT UPON THE SURFACE OF AN 

ELLIPSOID. 

By Thomas Craig, Fellow of the Johns Hopkins University. 

The ideas contained in the following brief paper were suggested to me 
in reading the twenty-eighth Chapter of Jacobi's Vorlesungen iiber Dynamik. 
A similar investigation may have been given before, but I have never seen 
anything on the subject. 

Let the point be acted upon by a force directed constantly towards the 

center of the ellipsoid and varying as the distance from the center. Let /?, a 

constant, denote the force at unit's distance from the center, and let a denote 

the force in the direction of the normal to the surface. We have then for the 

equations of motion of the point 

d 2 x x , a 

** = «■? +0* 



a) ?! = •-& + * 



*Z = a y~ 

dt 2 b 2 



d 2 z z , 



the ellipsoid being given by the equation 



a 2 ^ b 2 ^ c 2 ' 



Multiplying equations (1) by -^ , 2L 9 -^ respectively, and adding, we obtain 

CL C 

x d?x . y d 2 y z d 2 z _ [~ x 2 y 2 z 2 ~] 

^dt 2 ~ ] ~VdT 2 ' ] ~'?dT 2 ~^ a l~^ + '¥ + vi 



+45 + £+IJ- 



Calling Vp the reciprocal of the length of the perpendicular from the center 
on the tangent plane at the point x, y, z, we may write this equation 

, N x d 2 x , y d 2 y . z d 2 z , n 
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Differentiating the equation of the ellipsoid twice with respect to t, we have 

x dx . y dy z dz ~ 

a?~dl + b 2 dt + ~?di~ ' 
x d?x y d?y z_ cPz 

It 2 ~dt 2 + T 2 ~d7 z ^ r 7 2 dt 2 

1 (dx\ 2 1_ (dy\ 2 1 (dz\ 2 _~ 
+ ~a 2 \Jt) + ~¥ \dt) + ~? W " 5 
designate by P the last three terms of this last equation, then we have from (a) 

(2) * — P = op + (3 . 
Again, equations (1) give 

r 1 dx d 2 x 1 dy d 2 y 1 dz d 2 z~l V x dx y dy . z dz~\ 

\-a 2 ~dt~d¥ + l~ 2 dtdt 2 + 'd 2 di ~dt 2 \ ~ a Ltfdi + ¥di + V Tt\ 

B Vx_dx j^dy z_dzl m 

^ P la 2 dt^ b 2 dt^ c 2 dtJ ' 
this is simply 
/o\ dP dp 

(3) Tt= a dt' 

Eliminating a from (2) and (3), we have 

1. ^L _L- ^ ^ 

Jli + P + /3 dt ~ ' 
the integral of which is 

(4) f (P + (3) = A = const. 

Multiplying equations (1) by — - , -£ , and — respectively, and adding, gives 

< ^^,'^_^_i_dzd^z_ „r dx <ty dz~] 

dt dt 2 ^ dt dt 2 ~^ dt dt 2 ~ P L x di ~^ y Tt^~ * dl\ ' 
Integrating this we get 

^ ©V(£- + ©'=(£)'=^ + / + o + *. 

where B — const, of integration. For $ — or when no force acts towards 
the center, we have the simple case, 

s ■=. \/B . t + const. 
or the arc varies directly as the time. Any point on the ellipsoid can be 
given as the intersection of this surface with the two confocals 

* 2 + f I * 2 -i 
a 2 + a* ^ b 2 + Xx ^ c 2 + ^ ~ 

_f!_ + _i!_ + _^!_ _ i 

« 2 + A, ^ 5 2 + a 2 ^ C 2 + A, ' 
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Then we have X t and X 2 as the elliptic coordinates of the point, Xj_ = const, 
being the equation of one set of lines of curvature and X 2 = const, the equa- 
tion of the other set, the intersection of any two of which determines the 
position of a point upon the surface. We have now (Salmon's Geom. of Three 
Dimen., page 124), 

x* + f ■+ * = <* + *' + c 2 -(^ + a,) ; 

also, 

(6) tf s 2 =^ 2 +^ 2 +^=*^^ 

and for the differentials dx, dy, dz the known values 

2 ■•(<;» _«») («" — »») L V + »i V + xJ' 

and finally for the perpendicular from the center to. the tangent plane 

(8) * ^-+t + ^ = ^b =p . 

y } a 4 ^ b 4 ^ c 4 ale F 

From the expressions for dx, dy, and dz we can readily obtain. 

do? dtf dtf _ V-^ r <fr& «&§ 1 

W a 2 f J 2 + c 2 ~ 4 L(«2 +Xi )( & 2 + ^ iX<? 2 + ^ ) (a'+^Xft'+^X^+XoJ ' 

the first member of which equation, divided by dt 2 , is the quantity P, therefore 

nm ^ — a 2 Wf/ Vtf*/ 

* ; 4 L $ * 

where for convenience we write 

<D = (« 2 + ^)(& 2 + ^)(c 2 + ^i) 
<P = (a 2 + Xa) (b 2 + **) (<? + a,) . 

Substituting for P its value as obtained from equation (4) , this last becomes 

^ — a 2 fl /<ZXA 2 1 /^A 2 -| _ I? « 
t 11 ' ~~ T~ L^ VdtJ ~ ¥ V «fc/ J ~ ^ 2 _P 

where I) = obcA. Equation (5) becomes in elliptic coordinates 

(is) h=£± 1% (%)'-| (§)*] = c-p <*, + v) , 

where for brevity we write = /? (a 2 + 5 2 + c 2 ) + i?. Eliminate £ from 
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equations (11) and (12), and the result will be the differential equation, in 
elliptic coordinates, of the path of the point upon the ellipsoid : this is 

or 

/io\ V/l^tf/li , V A 2 «/l$ si 

in which the variables are separated. If again we make /? = 0, that is, if 
after the first impulse, no force acts upon the point but that in the direction 

of the normal, this equation becomes, on making — — = 2 , 

c 

= 0; 



V {a 2 + a*) (£ 2 + \) (c 2 -f a*) (^ + Ai) V (a 2 + a,) (* 2 + a 2 ) (c 2 + a 2 ) (0 2 + X,) 
this is the differential equation of a geodesic upon the ellipsoid. 

The integrals of these expressions will be elliptic for 6 = 0. If we change 
the variables \ and a 2 into new ones defined by the equations 

A 2 — — 

n 
our differential equation becomes 

(14) ^ 



V(a 2 £ + 1)(6 2 £ + l)(c 2 £ + l)(Df - ■<* + 0) 

J, ^ 

V (a 2 >7 + 1) (ftSz + 1) (<ftj + 1) (A 2 - 07 + 0) 

Write - = £and-^ = w, and call — £ l5 — £ 2 , — »h — >? 2 the roots of the 

quadratic equations £ 2 -f- #£ -f- m = , 

>7 2 + #»? + m — ; 

and further — - = «', — = b', — — c', and our equation can be written 



(15) 



S 1 



v(« + <o« + »o« + <o« + «,)« + &) 



/- * 



dh J rr const. 

V(>? + «')(>? + *')(>? + (»? + >7i)(>7 + *? 2 ) 



These are not in general reducible to elliptic integrals. If now we find a 2 as 
a function of a a , say x 2 —f{\) = /and *P = F(X X ) == F, then substituting in 
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either of the equations (11) and (12) we will be able to find the time if as a 
function of /^ . 

Making this substitution in equation (12), we readily find 



= r /[> 

d7n. 






When we can determine the values of f and F from equation (15), we 
will be able from this last to express the elliptic coordinates \ and 2, 2 as func- 
tions of the time, a problem which, for the case of a simple pendulum or 
point constrained to move upon the surface of a sphere, has been solved by 
M. Hermite in Orelle's Journal, vol. 85. 

We have seen that in general the path traced out by the moving point is 
not a geodesic, but was such a line where there was no central attracting force, 
or where /3 = 0. It is interesting, however, to observe that, for another value 
of /?, the point will move along a geodesic. Suppose that we have two quan- 
tities, t and a, the former either constant or a function of the time, the latter 
constant or a function of the arc s. And now make 

*'=[>+•!]£•= As 
* = [•+•£]&=** 

n r . ds~\ dz dz 

Our equations of motion now assume the forms 

d 2 x ax n dx 

W~a T + Pl Ji 

dt 2 ~ ' V "*" Hl 'dt 
d 2 z az , r> dz 

iw.-w^ ** x dt ' 

We have as before 

[" 1 dx d 2 x ' 1_ dy <Py_ 1 dz dh~\ _ [~f_ ^r # % _i_ £. §£\ 

La 1 did!* + ~¥ di dt 2 + ~d* dt dt 2 J ~ a U 4 dt V dt c 4 dU 

+A[i(|)VMI)VMI)l. 



or simply 



£=.£* + **. 
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And again, 

x d 2 x y d 2 y z d 2 z ___ fx 2 .y 2 .^ 






a 2 dt 2 ^ b 2 dt 2 ^ c 2 dt 



e+*+a 



r, fx dx . y dy , z dz~\ 

+ Pl La 2 di + Vdi + ?diJ' 

or simply 

— P = aj? . 
Eliminating a from these equations 

1 dp 1 tfP _ 9/? 

^W + P^zf-^ 1 ' 
from which we have 

*£ + ** = 20^=2 (tdt+o^-dt) . 
V P ^ dt s 

Wow writing 

2fvdtz= T 

2 Jotfs = 2 
we have as the integral of this equation 

We can also readily obtain from the equations of motion the following 

di 

-— — C"e }i< -' v + S) . 
dt 

These equations give 

__1_ /ds\ 2 _ pP 

C" 2 \dt) ~ C ' 

which becomes in elliptic coordinates 

o_ rv-^ pa (d^\ 2 _ ^ /d^\ 2 -j -i _ ^ aq-a 2 r_i at^a 2 _ i^ (d\\ 2 i 

6" 2 L 4 ~ U> \ <fr y qr V ^y J J — afo . 4 L<|> \dt ) *P V ^ / J ' 

abcC' 
Putting — — = — 2 , and multiplying through by dt 2 , we have after simple 

reductions 

v^r^i , VA4&4 



± ,/ 2 * =0, 



\/<& (a x + 2 ) V* (Xa + 2 ) 
or 

V (« 2 + ^) (& 2 + a,) (c 2 + a, x ) (a 2 + a, j V («» + a 2 ) (P + a 2 ) (<? 2 + a 2 ) (0 2 + a 2 ) ~ " 

This is the differential equation of the path of the moving point which is, 
as we see, a geodesic upon the ellipsoid. 



